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ABSTRACT 

Star formation happens in a clustered way which is why the star cluster population of a particu- 
lar galaxy is closely related to the star formation history of this galaxy. From the probabilistic 
nature of a mass function follows that the mass of the most-massive cluster of a complete 
population. A/max, has a distribution with the total mass of the population as a parameter. The 
total mass of the population is connected to the star formation rate (SFR) by the length of a 
formation epoch. 

Since due to evolutionary effects only massive star clusters are observable up to high ages 
it is convenient to use this 7\/,„ax(SFR) relation for the reconstruction of a star formation his- 
tory. The age-distribution of the most-massive clusters can therefore be used to constrain the 
star formation history of a galaxy. The method, including an assessment of the inherent uncer- 
tainties, is introduced with this contribution, while following papers will apply this method to 
a number of galaxies. 

Key words: galaxies: formation - galaxies: evolution - galaxies: star clusters - galaxies: 
stellar content 



1 INTRODUCTION 

During the last few years it has been recog nised that most an d 
. probably all stars form in embedde d clusters jLada & Ladall2003h . 
' The observational work, notably by'Larsen' 12002^, has established 
that star clusters ranging from the oldest globular clusters to the 
youngest low-mass objects have to be regarded as a continuous dis- 
tribution by mass — globular clusters are not fundamentally differ- 
ent from open clusters but merely the upp er mass end of the distri- 
bution. This has been shown explicitly bv lKroupa & Boilvl j2002l) 
on the example of the Milky Way population II spheroid. Today star 
clusters must be viewed as the "fundamental building blocks" of 
galaxies because they also determine the morphological appearance 
of wh ole galaxies through the physics of their formation l lKroupal 
l2005h . 

Star formation is therefore closely connected to the star cluster 
distribution in a galaxy and thus it is expected that the star forma- 
tion history of a galaxy leaves its imprint on the star cluster distri- 
bution. Since star clusters can live for a long time the star cluster 
distribution of a galaxy can be compared to a diary: Events like in- 
teractions of galaxies lead to an enhanced star formation activity. 
From this results a larger number of clusters being formed during 
the interaction time. This qualitative statement is well known and 
in this work we present a quantitative method to derive the star 
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formation history (SFH) of a galaxy directly from its star cluster 
content. 

Until now detailed SFHs can be determined only for galaxies 
which are at a distance allowing individual stars to be resolved, i.e. 
within the local group up to ~ 1 Mpc away. The colour-magnitude 
diagram (CMD) which is thusly obtained then allows the construc- 
tion of a SFH using theoretical isochrones. 

In distant galaxies, however, star clusters appear as compact 
sources whereas the individual stars give a homogeneous distribu- 
tion over the area. Modem instruments like the Hubble space tele- 
scope make it possible to obtain a cluster age- and mass-distribution 
for galaxies beyond th e distance where in dividual stars can be re- 
solved, e.g. for M SI jSastian et alj|2005|: cluster formation rate: 
iGieles et al.ll2005h and MIOI d Bianchi et a l."2005^. These galax- 
ies lie at a distance of ab out 7 Mpc (M51: Takats & Vinko 20d3; 
M 101 : lKeIsonetalJI 19961) which demonstrates the potential of our 
new method. 

A first approach to derive the star formation history from a 
cluster age- and mass-distribution could simply be to use all ob- 
served clusters in given time intervals. But this would lead to 
wrong results, clusters evolve and thus there are fewer clusters at 
higher ages. However, massive clusters evolve slowly, clusters with 
masses > 10^ have life ti mes comparable to a Hubb l e time 
teaumgardt & Makinol l2003h . IWeidner. Kroupa & Larserj ( l2004h 
found a relation which establishes the connection between mas- 
sive clusters and the galaxy-wide star formation rate: During a for- 
mation epoch a complete population of star clusters is formed, for 
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which the mass of the most massive cluster. A/max, depends on the 
star formation rate (the Mmax(SFR) relation). Consequently, the 
star formation history of a galaxy can be regarded as a sequence of 
such formation epochs, in which the most massive cluster carries 
the inform ation about the star fo rmation rate. 

Since IWeidner et al] j2004) used a deterministic law for the 
Afmax(SFR) relation, we briefly re-analyse their argumentation, 
allowing a distribution of the most massive cluster rather than a 
fixed value for a given SFR. Then we present the new method and 
conclude with testing it for some typical cases. 



2 COMPLETE POPULATIONS AND THE 
A/max (SFR)-RELATION 

2.1 Complete Populations 

The distribution of star clusters in a star cluster population is char- 
acterised by th e shape of the distribution function (shown to be a 
power law, e.g. IWeidner & Kroup j (12006) and references therein), 
the mass of the population, Aftot, and the mass limits, A/min and 
A^MAx. The cluster mass function can be written as 



= kM' 



A^^ 



where the normalisation constant k is determined by A^tot , 
and Mmax- 

In this context complete population denotes a statistically 
meaningful representation of the embedded cluster mass function. 
The underlying distribution function is defined on a mass interval 
appropriate to the mass of the population. This is a very general 
concept which is applicable to other cases where objects obey a 
distribution function, e.g. stars in a star cluster. 

The lower mass limit, A/min, is given by the physical mini- 
mum mass of a star cluster and is independent of the mass of the 
population. Until now, the physical minimal mass of a star cluster 
is not well known, here it is assumed to be 5 M© corresponding 
to grou ps of about a dozen stars such as those forming in Taurus- 
Auriga tericeno et al.ll200^ . 

For the upper mass limit. A/max, two cases have to be dis- 
tinguished: If a complete population has a mass which is larger 
than the limiting physical maximum cluster mass, then the distri- 
bution function is defined up to the physical maximum mass, a hith- 
erto not well understood quantity. Star clusters with masses larger 
than about IO^'Mq show complex stellar populations and are prob- 
ably distinct from th e "normal" star cluster content of a galaxy 
dWeidner et al] |2004) . Such arguments based on the structural 
properties of cluster s would imp l y the same physical maximum 
mass in all galaxies. iGieles et alj j2006h derived maximal cluster 
masses from the cluster luminosity function, obtaining masses be- 
tween 0.5 - 2.5 X lO'^Mo for NGC 6946, MSI and the "Anten- 
nae" (NGC 4038/39). Thus, while arguments exist for a limiting 
maximum star-cluster mass below which stellar populations are 
simp le (mon o-met allic a nd -age), Mieske, Hilker & Infante ( 2002) 
and [Martini & Hoi ( |2004|) show that ultra-compact dwarf galaxies 
(Af > lO^M©) may be an extension of the "cluster" formation 
process to large masses. Consequently, we do not limit the "clus- 
ter" masses but allow these to formally reach A/max ~ lO^M© 
for sufficiently high SFRs. 

The case that Aftot is smaller than the physical maximal mass 
implies that no star clusters more massive than Aftot can exist in 
this population. Therefore Aftot is the upper mass limit of the clus- 
ter mass function. This includes the case that a population can con- 
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Figure 1. M asses of the b rightest clusters vs. presen-day SFR in galax- 
ies, data from L^^eiJ j2002h (absolute magnitudes converted to masses by 



(1) IWeidner et"ai] 



2004) . The lines shown with the data ai'e based on our sta- 



tistical point of view, presented in Sec. 12.21 The parameters St and /3 ai'e 
chosen to fit to the data with St = lOMyr, p = 2.4, and A/min = 5M0. 
The solid line is the expectation value for the distribution of the most mas- 
sive cluster (eq. |6) in dependence of the total mass of the population. The 
dashed lines are the borders of the region in which 2/3rd of all most massive 
clusters are expected, calculated using eq.|5] 



sist only of one cluster with the mass A/tot, but this is very improb- 
able. 

Thus by using the total mass as the normalisation criterion. 



Mi^i{M)dM, 



a/min 

the normalisation constant becomes 

A/tot(2-/3) 



k = 



A/, 



2-/3 



A/, 



2-/3 ■ 



(2) 



(3) 



-"MAX ^"MIN 

The total number of clusters in a population, A'tot, follows from 
iVtot (A/tot) = [^''t.t^ \ ci{M)dM. 

IWeidner et al.fj2004l) argued that a complete population of star 
clusters is not made up by all clusters ever formed in a galaxy, but 
by a subset of clusters formed during a. formation epoch. Assum- 
ing that all stars form in star clusters the total mass of a complete 
population is then given by the product of the SFR and the length 
of the formation epoch. 



Aftot = SFR X 5t. 



(4) 



Thus, given a certain (short) 5t, the total mass and thereby implic- 
itly the upper mass limit of the distribution function depend on the 
current SFR. With this description massive clusters can only form 
if there is much star forming activity, while quiescent phases only 
produce low-mass clusters. 

2.2 Tiie Afmax(SFR) relation 

Observations give evidence that the brightness of the brightest clus- 
ter in a galax y depends on the present-day SFR in the galaxy 
( lLarsenll2003) . Since the brightest cluster in a galaxy is usually 
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Figure 2. The distribution of Mmax for a population with Mtot = 
10*Mq, P = 2.4, and AfjvilN = SMq. The number of A/max expected 
to lie below and above the average A/max is different due to the asymme- 
try. The median, A/i/2> a different value than the average. 2/3rd of all 
Afmax are expected in the region delimited by lines marked with 1/6 —> 
and^ 1/6. 



young it can be inter preted as the most ma ssive cluster of the cur- 
rent formation epoch. Tweidner et afl ^2004) converted the luminosi- 
ties to masses, as shown in Fig. 1. The data show a large scatter 
which in our interpretation results from the distribution of the most 
massive cluster. Since the distribution of A/max for a given A/tot is 
known the data can be used to determine the length of the formation 
epoch. 

In our description of star cluster populations the upper 
mass limit, A/max and the most massive cluster, Afmax, are not 
identical. For an ensemble of populations with the same total 
mass, Mmax has a distribution parametrised by Mtot and fj (cf. 
lOev & Clarkel2005l) . This distribution can be written as a probabil- 
ity density, 



<^(Mmax) 



1 

Wo 



U{M)dM 



JVto 



Ccl(Mma0) 



where A'^tot and the nor malisation of ^^ i depen d on Mtot . This is 
different to the ansatz of IWeidner et alj ( |2004|) where A/max was 
assumed to be identical for all populations with the same A/tot, i.e. 
not distributed. The distribution of Mmax is asymmetric because of 
the asymmetric cluster mass function and is characterised by the 
average mass of the most massive cluster. A/max, given by 

-^-^MAX 

(6) 



M 



A'/max ( A^max ) dM[-,-, 

ax ■ 

Figure|2]shows the distribution of Mmax and the location of A/max. 
Due to the asymmetry the median. Mi/ 2. does not have the same 
location as the average, Afmax, but lies below it. Therefore it is 
expected that in an ensemble more Afmax lie below A/max than 
above. 

The observations shown in Fig. [T] are an ensemble of most 
massive clusters of populations with different total masses. Since 
<^(Mmax) is parametrised by Mtot it is possible to derive the av- 
erage A/max in dependence of A/tot and via eq. |4] also in depen- 



dence of the SFR. This analytical curve could be compared to a line 
derived from observations e.g. by least squares fitting. The least 
squares fitting procedure would give a line that leads through the 
region where most of the A/max are, but this does not match A/max 
due to the asymmetry of </f>(Mmax). Most of the Afmax lie below 

A/m ax. 

To constrain 5t from the observations it is more convenient to 
use the region where a certain fraction of the data is expected. The 
location of this region - at high or low values of Mmax - depends 
on the duration of the formation epoch: A long formation epoch 
pushes the region towards high masses, a short formation epoch 
to lower. A larger /? steepens the relation. The best fitting region 
where 2/3rd of the m ost massive clusters a re expected is shown 
in Fig.[T] As found bv lWeidner et al.l j2004h the values of St = 10 
Myr and /3 = 2.4 provide a good fit to the data. Complete popula- 
tions of star clusters form in formation epochs lasting for 10 Myr 
which is comparable to the time -scale of th e emer gence of clus- 
ter populations from spiral arms tegusa, Sofue & Nakanishill2004l: 
iBonnell et ai]|2006h . 



3 A METHOD TO DERIVE STAR FORMATION 
HISTORIES USING STAR CLUSTERS 

3.1 Concept and restrictions 

The analysis of the Afmax(SFR) relation suggests that the SFH 
of a galaxy can be interpreted as a sequence of formation epochs 
lasting for ~ 10 Myr. A star cluster population with a total mass de- 
termined by the current SFR emerges during each formation epoch. 
Therefore it should be possible to infer the SFR from the properties 
of this population, but for the largest part of the lifetime of a galaxy 
only a small fraction of coeval clusters is observable making the 
determination of the population mass difficult. 

On this account we propose a different approach using only 
the most massive clusters of each formation epoch. The mass of the 
most massive cluster also depends on the current SFR, but the dis- 
tribution of them does not allow the derivation of the SFR from only 
one most massive cluster. If it is assumed that the galaxy-wide SFR 
changes significantly only on a time-scale that includes a number of 
formation epochs, then the clusters of this set of formation epochs 
can be seen as an ensemble of identical cluster populations. For 
the ensemble average of the most massive clusters the mass of the 
population can be calculated using the Mmax (SFR) relation: The 
probability density (j!)(A/max) is parametrised by Aftot = SFR x 5t 
and [3 and consequently also A/max, written symbolically as 



Mmax = /(SFR), (7) 

where /(SFR) is given by the integral of eq.|6] Thus the SFR cor- 
responding to this ensemble follows by inversion of the previous 
equation, 

(Afmax). 



SFR 



/" 



(8) 



The inversion of the integral (/) is done numerically. 

This gives the general idea to reconstruct SFHs: The lifetime 
of a galaxy is divided into time windows containing a number of 
formation epochs. For each of the windows the average mass of the 
most massive clusters is calculated and from this mass the underly- 
ing SFR is derived. 

The age determinations of ava ilable cluster data for gal axies, 
as e.g. the Large Magellanic Cloud jde Griis & Anderil2006h . usu- 
ally have uncertainties of ~ 0.4 dex. If the observational situation 
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is optimal, i.e. observations in the most suitable fi lters could be 
made, then the age uncertainty can be much smaller. Ide Griis et al.l 
1 I2OO5 I) obtained A log(age/yr) < 0.15 ' 'in the majority of cases" 
for conditions as for NGC 3310. Since it does not make any sense 
to try to detect variations of the SFR on time-scales shorter than 
the age uncertainties of the data, the averaging window has to be 
chosen to have the same length (or longer) than the available age 
uncertainty. Because the error is constant in logarithm the length 
of the averaging window depends on the age. We chose a length 
of 0.5 dex for the averaging window. A sequence of neighbouring, 
independent averaging windows leads to a SFH determined at dis- 
crete points in time over the lifetime of a galaxy. For a continuous 
SFH the averaging window is moved in 10 Myr-steps. At old ages 
the length of the averaging window is reduced to ensure that the 
oldest formation epoch used in the averaging window contains a 
cluster. An averaging window filled only halfway would lead to a 
systematically underestimated SFR. 

As the averaging window increases with age, the minimum 
duration of resolvable events in the SFH also increases. This is an 
inherent restriction of our proposed method. Short bursts of star 
formation which happened at large ages cannot be resolved. 



3.2 Corrections for dynamical cluster evolution 

The cluster ages and masses for large samples of clusters are usu- 
ally determined by fitting models to observed spectral energy distri- 
butions. The age and the "initial" mass of a cluster are parameters 
for the fitting routine. In all cases the cluster mass is determined by 
scaling a model "initial" mass. Since the applied models usually 
only consider mass loss due to stellar evolution and not due to dy- 
namical evolution, the fitted mass of a cluster does not correspond 
to the initial cluster mass for which the A/max (SFR) relation is 
valid. 

Star clusters are no isolated static objects: the gravitational 
force keeps the stars in constant motion with respect to the centre of 
mass, which itself moves on an orbit through the host galaxy. Stars 
can evaporate from the region dominated by the cluster potential 
and leave their star cluster. This mass loss of the cluster due to 
dynamical evolution depends on the eccentricity of the orbit and 
the distance to the galactic centre. Furthermore, it depends on the 
mass of the host galaxy which determines the strength of the tidal 
field. 

Since these parameters are mostly not available, the analytic 
model from lLamers et al. I( l2005l) for the averag e mass loss suffered 
by a cluster in a galaxy is used. iLamers e t al. (2005) find a good 
agreem ent of their model with the results of Baumgardt & Makinol 
( l2003l) . who give a formula derived from A*'-body experiments. 

For the statistical model adopted here it is assumed that cluster 
disru ption depends only on the initial mass jBoutloukos & LamersI 
I2OO3 ). In this case, the initial mass can be calculated from the ob- 
served mass and age using 

where 7 = 0.62 is identical for all galaxies, to describes the tidal 
field and can be determined from the dissolution time of a IO^Mq 
cluster, 

Veeoy 

The parameter ti has been determined for a number of galaxies (cf. 



Boutloukos & LamersI l2003l : iLamers. Gieles & Portegies ZwartI 
2005b . 

3.3 Upper and lower limit for the SFH 

The observed cluster content of a galaxy usually does not provide 
a cluster for every formation epoch. The number of observed clus- 
ters older than a few Gyr is much s maller than the number of for- 
mation epochs. Ide Griis & Anderj ( l2006h found only « 10 clus- 
ters older than 4 Gy r in the Large Magellanic Cloud, as similarly 
iBastian et alj ( l2005h for M5I. This originates from a SFR which 
was so low that no clusters were formed being massive enough to 
be visible today. The brightness limit of the observations is there- 
fore an upper limit for the brightness of the actually formed most 
massive cluster. With a cluster evolution model the brightness limit 
at a given time can be converted to an initial mass which then is the 
upper mass limit fo r the most massive cluster. The GALEV models 
jSchulz et alj|2002l) give the luminosity evolution J^{t) of a "sim- 
ple stellar population" (i.e. single burst, single metallicity) with a 
mass of 1.6 x 1O^M0 for different photometric bands. By scaling 
M{t) to the limiting magnitude of the observation, A^nm, with the 
band chosen according to the bands u sed in the cluster-m ass fitting, 
the limiting mass can be derived (cf. iHunter et al.ll2003h : 

Mita(t) = 1.6 X 10^+°-^(^W-^"»)Mq. (II) 

The GALEV models do not take dynamical evolution into account, 
therefore A/nm has also to be corrected for dynamical evolution as 
described in the previous section {M = in eq.|9). 

Now the SFH can be derived with as the most massive 
cluster in those formation epochs that do not contain any clusters, 
leading to the upper estimate of the SFH. Since A/nm increases for 
older ages due to the internal evolution of clusters, the derived SFR 
also increases, which does not necessarily reflect the underlying 
SFH. 

The lower limit of the SFH is given by using Afmax = in 
the epochs containing no cluster, corresponding to the assumption 
that during these epochs no star formation took place at all. 

3.4 Self-consistency checks 

In the next section tests to analyse the systematic errors of our 
method are made using modelled SFHs. For an observed nearby 
galaxy the obtained results can be compared to independently de- 
termined quantities. The SFH derived using star clusters should be 
similar to the one obtained using colour-magnitude diagrams. This 
comparison can be difficult if the used regions of the galaxy differ. 

The total mass of stars in a galaxy can be calculated since 
with our method a SFH for an entire galaxy is derived. For each 
10-Myr epoch a SFR has been determined, from which the mass of 
the formed stars can be calculated. The sum over all epochs gives 
the mass of the stellar content of the galaxy. When deriving this 
mass, stellar evolution is not taken into account, i.e. all stars that 
ever formed are counted regardless of whether they still exist or 
not. This can be compared to independent determinations of the 
stellar content of a galaxy. 

Furthermore, the cluster formation rate should reflect the 
structures found in the SFH. For each most massive cluster an ap- 
propriate total number of clusters should exist. This comparison of 
the cluster formation rate (i.e. number of clusters per time) and the 
SFH can be done e.g. by using a cluster formation rate also derived 
using a moving time window. 



Star formation histories of galaxies 5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I ■! I I I I ^ I I I I I 




" i I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



Pi 




Pi 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



• V 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



_2 I I I I I I I I I 

7 8 9 10 

logio Age/yr 

Figure 3. Test of the new method for a modelled constant input SFR 
(dashed line in the lower part). In the upper part the most massive clus- 
ters of each 10-Myr formation epoch are shown (dots). The solid line in the 
lower part is the reconstructed SFH using the method described in the text. 

4 MODELLED STAR FORMATION HISTORIES 

To test our method synthetic star cluster populations were gener- 
ated for different SFHs. The aim is to verify if an input SFH can 
be re-extracted and to study the effects of the averaging. For this 
purpose we first consider the simplest case with optimal data, i.e. a 
constant SFR and no measurement uncertainties for the age. Mod- 
els with a varying SFR that include the age uncertainties show the 
capacity of our method. For clarity and to focus on the effects of 
the averaging and age uncertainties, cluster evolution is only con- 
sidered in the last model (Fig.|9j. 

The general procedure of our models is as follows: In time 
steps of 10 Myr complete cluster populations are generated with 
a total mass given by eq. |4] and the mass limits Mmin = 5M0 
and Mmax = min(Aftot, MMAx,phys) with MMAx.phys = 
IO^Mq. Then an age uncertainty is assigned to each cluster 
with age r, drawn randomly from a Gaussian in logarithmic age, 
A/'(logj^Q r'; ^ = logifl r, Or)- The values chosen for the variance 
Gt are 0.15 dex and 0.35 dex, correspo nding to the typical u ncer- 
tainty range of the SED fitting method I de Griis et alj|200^ . Un- 
certainties larger than 0.5 were rejected and generated again until 
they are smaller than 0.5. The procedure for synthetic cluster evo- 
lution is described in the Section of the respective model. For each 
10 Myr interval only the most massive cluster is shown in Figs.[3j- 

s 



4.1 Constant SFR 

Figure [3] shows the distribution of the most massive clusters for a 
constant input SFR (dashed line in the lower part) not including 
age uncertainties. It appears as if the masses of the most massive 
cluster increase with age, which is a size-of-sample effect. Due 
to the logarithmic axis 90 % of the clusters are in the age range 
9 ^ logj^o ^ 10' therefore the probability to get a very massive 
cluster by chance is higher in this range. Similarly the probabil- 
ity to sample clusters in the mass range below the average mass of 
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Figure 4. Same as Fig. [3]but for a slowly linearly decreasing input SFH. 
Here an age uncertainty with cr-r = 0.15 dox was included, leading to the 
decline at high ages as discussed in the text. 

the most massive cluster (Mmax, eq.|6j increases. This leads to the 
wedge-like shape of the distribution of the most massive clusters. 

The derived SFH (solid line) only shows a large variation at 
young ages, where the large averaging window of 0.5 dex in log 
age contains only a few 10 Myr formation epochs. Therefore the 
reconstructed SFH is calculated from too few clusters. For ages 
older than ~ 100 Myr the reconstructed SFH agrees within the 
expected statistical variations (as discussed in Section |431 l with the 
input SFH. The features in the reconstructed SFH at an age of 10 
Gyr are artefacts due to the shrinking averaging window. 

4.2 Varying SFR 

Figure |4] shows the SFH obtained from a linearly decreasing SFR. 
In this model an age uncertainty with Or ~ 0.15dex was in- 
cluded as described above. This causes a systematic deviation of 
the reconstructed SFH towards smaller SFRs at old ages. Although 
cluster populations were only generated up to an age of 10 Gyr 
the ages of individual clusters can be allocated up to maximally 
logj^Q r = 10 + 0.5. Thus, for many old formation epochs not 
the actual most massive cluster corresponding to it is used but the 
second or third etc. most massive. Therefore the SFR is underesti- 
mated. Besides this effect introduced by the age uncertainties the 
behaviour of the reconstructed SFH is similar to the one of the pre- 
vious model: The reconstructed SFH follows the slow change of 
the input SFH within the same degree of deviations. There is a large 
scatter in the SFH for ages younger than ~ 100 Myr. Also as in the 
constant case the artefacts at the oldest ages (logj^Q r > 10.3) are 
visible. This model allows to conclude that our method is capable 
of reproducing slowly changing SFHs. 

The modelled SFH shown in Fig. |5] has three well-separated 
bursts lasting from 0-200 Myr, 1 000-1 200 Myr, and 8 800-10 
000 Myr. During the bursts the SFR is increased by a factor of 10. 
Again the cluster ages were generated with an age uncertainty with 
(Tr = 0.15 dex. Our method results in an undulating SFH with 
peaks roughly coinciding with the centres of the bursts. The oldest 
burst is affected by the decline of the reconstructed SFH due to the 
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Figure 5. Same as Fig.|4]but for a SFH with tliree well-separated bursts. 

way how the age uncertainties are assigned. Thus the maximum 
is shifted towards younger ages. In the reconstruction the shape 
of the SFH is much less pronounced than in the original. The age 
uncertainties and the averaging procedure lead to a larger width of 
the recovered bursts. As the most massive clusters of a burst are 
spread over an interval longer than the burst the amplitude of a 
recovered burst also decreases. 

In Figures[3]to[5]the input SFHs guide the eye to see the appro- 
priate structure in the SFHs derived from the cluster distributions. 
The distribution of the most massive clusters, the age uncertainties, 
and the averaging process in our method lead to a much less dis- 
tinct structure in the reconstructed SFH than in the initial one. In 
Section l4!4l we present a criterion to decide which features in a re- 
constructed SFH are caused by variations in the initial SFH and are 
not merely due to systematic effects. For this purpose it is neces- 
sary first to investigate the systematic effects of our method, which 
is the object of the next Section. 



4.3 Statistical scatter in tlie reconstructed SFH 

The models presented above show that the input and the recon- 
structed SFHs differ, especially for ages younger than ~ 100 Myr. 
Due to the small number of clusters used for the reconstructed SFH 
in the respective formation epochs the scatter in the reconstructed 
SFH increases. To investigate the expected scatter in the recon- 
structed SFHs, a sample of star cluster populations with the same 
SFH was created. For each of the input SFHs, 1000 synthetic clus- 
ter populations (i.e. 1000 galaxies) were generated with different 
initial random seeds and the SFH was re-extracted from the clus- 
ters. Then the sample average, SFR(t), was calculated for each 
formation epoch, 

1000 

1=1 

where SFRi (t) is the SFH of an individual cluster population i. To 
achieve an estimate of the statistical spread we calculated the pos- 
itive and nevative deviation, (Tsfr+ and asFR-, from the average 
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Figure 6. Averaged reconstructed SFHs (thick solid lines) of a sample 
of 1000 synthetic cluster populations for different input SFHs (thin soHd 
lines). Also shown are the corresponding deviations, o"sfr+ and (Tsfr_, 
as discussed in the text (dotted lines, eq. ll3l and ll4t . 

a) constant input SFH, no age error. 

b) three well separated bursts, age error cft = 0.15 dex. 

c) three well separated bursts, age error cft = 0.35 dex. 

d) one burst per Gyr, age error a-r = 0.15 dex. 



SFR at a given time, 

'^sFR+(t) = |SFR.(t)-SFR(0|3) 

SFRi(t)-SFR{t)>0 

asFR-W = |SFR,(t)-SFR(a^) 

SFRi(t)-SFR(t)<0 

where N± [t) is the number of SFRs larger or smaller than SFR(t) 
at a given time t. This particular choice of an individual positive 
and negative average deviation will be discussed after the average 
SFHs. 

Figure [6] shows the results of the experiments. For the experi- 
ment presented in panel a) a constant SFH without age errors was 
chosen. The sample average and the input SFH agree well, differ- 



Star formation histories of galaxies 7 



Table 1. Original and reconstructed total masses of the stellar content for 
the different input SFHs. Output A/tot is the ensemble average of the 1000 
galaxies with the average deviation. 



10 



Model: 



input Aftot 
[logio Mq] 



output Mtot 
[loglO M©] 



Fig.|2]a) Constant SFR 
no age error 

O.OlM0/yr 8.00 8.02 

O-lMg/yr 9.00 9.01 

IMg/yr 10.00 10.01 

lOMg/yr 11.00 11.01 
Fig.|7]b) Three bursts 

ar = 0.15 dex 9.39 9.40 
Fig.|7]c) Three bursts 

(Tr = 0.35 dcx 9.39 9.41 
Fig.|2]d) One burst / Gyr 

ar = 0.15 dex 10.45 10.48 



+0.02 
-0.02 

+0.03 
-0.03 

+0.04 
-0.04 

+0.05 
-0.04 



+0.05 
-0.04 



+0.06 
-0.04 



-0.07 
-0.05 



ences occur only for young ages. This is because of the small num- 
ber of formation epochs used for reconstructing the SFH at these 
ages. Since the averaging window contains few epochs the scat- 
ter in the reconstructed SFH increases, which is also visible in the 
progress of (Tsfr±. The average deviations decrease with age be- 
cause the number of formation epochs used for the reconstructed 
SFH increases as a consequence of the averaging window which 
moves in logarithmic time. At the oldest ages the average deviation 
increases again since the averaging window contains a decreasing 
number of formation epochs with a cluster, which is the same effect 
as for young ages. As described in Section|4T|the SFH ends at the 
oldest ages with an artifact. 

Panels b) and c) of Fig. [6] show the results for the SFH with 
three bursts and different age uncertainties. For the smaller age un- 
certainty, (Jt = 0.15 dex, all three epochs of enhanced star forma- 
tion can be identified, i.e. the variation of SFR(t) is comparable or 
larger than (7sfr± {t). However, the reconstructed shape of the sec- 
ond and third burst is much wider and less pronounced. The larger 
age uncertainties (ctt = 0.35 dex) lead to a reconstructed SFH 
where the SFR decreases with age and only allow the reconstruc- 
tion of the youngest burst. Due to the age uncertainties SFR(t) 
drops at old ages, as discussed in Section |42l The average devia- 
tions crsFR± behave similarly to the constant case. 

In the case of one burst every Gyr, lasting for 200 Myr, (Fig.[6l 
panel d) with the small age uncertainty (ctt = 0.15 dex) only 
the youngest burst can be detected. Then the reconstructed SFH 
declines until a minimum between the first and second burst is 
reached. From the second burst on only a constant SFR with an 
intermediate value can be recovered. (7sfr± has the same features 
as before. 

A comparison of the averaged recovered SFHs and the input 
SFHs of these models show the capabilities of our method. Due to 
the averaging process structures in the SFH on time scales shorter 
than the averaging window cannot be reconstructed. In the other 
cases the sensitivity for structures in the SFH depends on the qual- 
ity of the age determination. For the situation of a SFH only slowly 
varying or with bursts that are well separated and small age errors, 
the shape of the underlying SFH can be extracted with our method. 
Large age errors of the star clusters or highly variable SFHs do not 
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Figure 7. Relative errors (thin solid lines) for different constant SFHs with 
a SFR of O.OlMo/yr, O.lMo/yr, IMo/yr and lOMo/yr and for the 
bursting cases as above, Fig.|6] The thick solid line shows an upper envelope 
of the relative error as described in the text (eqs. ll6l and ll7) . 



allow us to recover all of the initial features in the derived SFH. 
However, even in these cases the absolute value of the recovered 
SFR is of the same order of magnitude as the actual one. 

Thus, the total recovered stellar mass (A/tot, rcc = 
J SFR(t)dt) corresponds to the underlying SFH, summarised in 
Table [T] The ensemble average of Mtot, rcc (average of the 1000 
galaxies) equals to the input value and no bias is introduced. 

The second aim of the models is to investigate the scatter in 
the reconstructed SFHs. Two properties of the average deviations, 
crsFR+ and ctsfr-, are noticeable: First, in the log-log diagram 
the positive and negative deviation have for a certain age approxi- 
mately the same distance to SFR(t). A reason for this effect could 
lie in the asymmetric distribution of the most massive star cluster 
(Fig-dJ. Thus, using the average of o-sfr+ and ctsfr- would lead 
to a wrong uncertainty estimate. The second property of both aver- 
age deviations is that their values relative to SFR are independent 
of the value of SFR. This is visible in Fig.|7] where the quantities 



0"SFR±(i) = 



g"SFR±(^) 

SFR(t) ' 



(15) 



the relative average deviations, are plotted. To show the inde- 
pendence of (TSFRi from the shape of a SFH the SFHs from 
Fig. [6l panels b) - d), are used. The model with a constant 
SFR of O.lMo/yr (Fig. (6] panel a)), and additional models with 
O.OlM0/yr, lM0/yr, and IOM-q/jx demonstrate the invariance 
from the absolute value of the SFR. The independence of the rel- 
ative average deviations from the shape and absolute value of the 
SFH makes this quantity suitable for estimating the uncertainties 
of the method. For the implementation in our method for deriving 
SFHs of individual galaxies we used the analytic fitting formulae 



o"sfr+ (r) 
and 

o"sfr- (t ) 



45 



1 + expl.3(logio^ - 5.4) 



l + oxpl.2(logi(, r-8.0)' 



(16) 



(17) 
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Figure 8. Reconstructed SFHs from Figs.[3]|4] and|5](top to bottom: con- 
stant, linearly decreasing, with three bursts; thick lines), shown with the 
null hypothesis of a constant SFH (thin solid lines, (SFR)). For a signif- 
icant variation of the reconstructed SFH it has to cross the Icrgpf^j- line 
(inner dotted lines). The outer dotted lines are at 2(tsfr±. 



These fits are a conservative estimate and lie slightly above the ex- 
perimental data. In the method the absolute average deviation is 
then calculated by 



'7Sfr±(t) = (?sfr±(t) X SFR(r). 



(18) 



It will be used in the criterion to detect significant variations in a 
SFH, discussed next. 



4.4 A criterion to detect significant variations in a SFH 

Due to the probabilistic distribution of the most massive cluster 
scatter in the reconstructed SFH is expected. As visible in Figs. |3] 
and |4] the reconstructed SFH can mimic periods of reduced or en- 
hanced star formation. However, the variations generated in this 
way are not caused by real events. The average deviations derived 
in the previous Section can be used to give the region of SFHs 
compatible with the reconstructed one. But to verify variations in a 
SFH we suggest to disproof the hypothesis that there are no varia- 
tions. This is done by setting up the null hypothesis of a constant 
SFH with a SFR equal to the time-averaged SFR, (SFR), of the re- 
constructed SFH. If the reconstructed SFH leaves significantly the 
lo"SFR± -region of the null hypothesis, a constant SFH, can hardly 
be supported. Thus the actual SFH of the galaxy has to have varia- 
tions. This is our suggested criterion of significance. Note that this 
criterion implies that troughs and maxima that result from a truly 
variable SFR differ by more than 2ctsfh,±- 



(SFR) is given by the integral over the SFH, divided by the 
time. With our method two SFHs are reconstructed, the lower limit 
with gaps where no cluster was observed, and the upper limit with 
the fading limit mass used in the gaps. Therefore (SFR) has to be 
calculated as the mean of the average SFR of both limits. In reality 
these limits start to differ for older ages, caused by incomplete- 
ness due to cluster evolution and observational limits. Especially in 
galaxies with strong cluster evolution the upper limiting SFH can 
shift (SFR) towards unreasonable values if it is integrated over all 
ages. To prevent this we integrate only up to the age where the log- 
arithms of the upper and lower limit differ by less than 0.2. 

The outcome of this procedure is displayed in Figs. [S] and |9] 
Figure [8] shows again the SFHs of Figs. [3] El and|5] now with the 
criterion for variations. Since in these models no cluster evolution 
was incorporated, we integrated up to an age of 6 Gyr (where the 
thin solid line stops) to obtain (SFR). In the constant case the input 
SFR (O.lMe/yr) and the average ((SFR) = O.12M0/yr) are in 
good agreement. The extremes of recovered SFH barely exceeds 
the lo"sFR± region. Thus the null hypothesis of a constant SFH 
cannot be rejected, as is correct in this case. 

The reconstructed linearly decreasing SFH lies, during most 
times, far away from the average. At young ages the recovered SFH 
lies sometimes even below 2(7sfr,-, and before the decline due to 
the age uncertainties at old ages it rises above 2(jsfr+- From this 
can be deduced that the null hypothesis is not consistent with the 
reconstructed SFH. A clear trend of increase with age is visible. 

The bursting case is harder to identify. Only the first burst (0- 
200 Myr), and the last dip (1.2-8.8 Gyr) leave unambiguous traces, 
albeit the shape is less pronounced than that of the underlying SFH. 
The two other bursts and the first dip merely allow an "educated 
guess" of the shape of the actual SFH. As above a constant SFH 
is not consistent with the reconstructed SFH. However, the exact 
structure of this SFH cannot be determined with sufficient certainty, 
although it is visible. 

The above cases show that our suggested criterion of signifi- 
cance allows us to distinguish between features in the actual SFH 
and artefacts due to the method. 



4.5 Constant SFR including cluster evolution and an 
observational limit 

The models described above do not account for cluster evolution 
and the observational limiting magnitude for cluster detection. To 
show the consequences of these effects a galaxy-model with condi- 
tions similar to the Large Magellanic Cloud (LM C) was generated . 
The stellar mass of the LMC was determined bv lKimetal.Ul998l) 
to be 2.0 X 1O^M0. Assuming a constant SFH leads to a SFR of 
« O.IM0 /yr, which we used in our model. The cluster-disr uption 
parameter is = 7.9 x lO^yr jBoutloukos & Lamers"2003) and 
the flux limit is A4iim = — 3.5mag jHunter et al...2003) in the V 
Band. The generated data and the results are shown in Fig.[9l 

Because of the weak tidal field clusters evolve only slowly and 
the high mass-end of the clusters is similar to the case without clus- 
ter evolution. The effect of the observational limit is clearly visible 
as a cut-off in the lower part of the cluster distribution. The min- 
imum observable mass increases with time due to the luminosity 
evolution of the clusters. 

As in the previous models the reconstructed SFH deviates 
from the input value at young ages because of the small number of 
formation epochs used for averaging. Then there is a period of good 
agreement, until some formation epochs contain no clusters any 
more. Due to the cluster evolution the clusters of these epochs are 



Star formation histories of galaxies 9 



-1 



-2 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I' 




I I I I I I I kij I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r 




1 1 •' 



9 



10 



logio Age/yr 



Figure 9. Recovered SFH for a LMC-type cluster model that includes dy- 
namical evolution and the missing of clusters due to the observational flux 
limit (dashed line in the upper part). As previously the upper panel shows 
the cluster distribution, whereby here each cluster is shown with its current 
mass after evolving it dynamically, and the lower panel contains the recon- 
structed SFH (thick Hues). Due to missing clusters the reconstructed SFH 
branches in two parts, as described in Section l431 The thin soUd line is the 
constant SFH ((SFR,))used to detect significant variations of the true SFH 
and stops when the averaging is stopped. Dotted lines indicate the Icgpjj^^ 
and 2(TSFH-|- regions. 



dissolved or have lost such a large fraction of their stellar content 
that they cannot be detected. Therefore the reconstructed SFH now 
shows two branches corresponding to the upper and lower limit. 
Assuming that no detection of a cluster means that there was no 
star formation activity leads to the lower limit. Using the detection 
limit as the mass of the most massive cluster gives the upper limit. 
The true SFH lies between both limits, which is indeed confirmed: 
Until the branching of the upper and lower limit the 1(7sfr,± -region 
contains the reconstructed SFH. From this point on only a rough 
estimate of the SFH can be obtained. From the input SFH a total 
stellar mass of 1 X 10^ Mq was built up. The upper and lower limit 
reconstructed 0.78 x IO^Mq and 1.74 x lO^M©, embracing the 
model value. That the reconstructed mass range contains the known 
input value constitutes a consistency check. 



5 CONCLUSIONS AND SUMMARY 

Based on the assumption that all stars form in star clusters it is 
possible to explain the relation between the brightest clusters in a 
galaxy and the present star formation rate, assuming that the bright- 
est cluster is also the most massive. During a formation epoch last- 
ing for r; 10 Myr a complete population of star clusters is formed. 
The mass of the most massive cluster obeys a distribution func- 
tion that can account for the scatter in the A/max (SFR) relation. 
The SFH of a galaxy can then be seen as a sequence of formation 
epochs. Starting from this we presented a new method to derive 
SFHs using star clusters taking into account the statistical proper- 
ties of the most massive cluster as well as their dynamical evolu- 
tion. 

The method was tested for a number of model SFHs for which 



synthetic cluster populations were created. The tests show that our 
method is capable of reproducing the modelled SFHs if they are 
only slowly varying or have bursts which are well separated. To be 
resolved, the time between two short-time bursts needs to be longer 
than 0.5 dex, the time over which is averaged. However, the typical 
uncertainties in the age determination and the need for averaging 
do not allow a shorter averaging window. Artefacts result from av- 
eraging over too few small age bins and from missing data at high 
ages. The example SFHs show to which degree our method can be 
used to make confident statements about the SFH of a galaxy. 

A model including realistic conditions for observation and 
cluster evolution also leads to good agreement between the input 
and the reconstructed SFH. In following contributions we will ap- 
ply this method to the galaxies LMC, SMC, M51 and MlOl. For 
the LMC we will compare this new method to the results obtained 
using the CMD method. 
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